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Abstract 



A classical microscopic theory of rovibrational motion at high angular momenta in sym- 
metrical non-linear molecules AB2 is derived within the framework of small oscillations near 
the stationary states of a rotating molecule. The full-dimensional analysis including stretch- 
ing vibrations has confirmed the existence of the bifurcation predicted previously by means 
of the rigid-bender model [see B.I. Zhilinskii and I.M. Pavlichenkov, Opt. Spectrosk. (USSR) 
64, 413-414 (1988)]. The formation of fourfold energy clusters resulting from the bifurcation 
has been experimentally verified for H2Se and it has been demonstrated in fully-dimensional 
quantum mechanical calculations carried out with the MORBID computer program. 

We show in the present work that apart from the level clustering, the bifurcation produces 
physically important effects including molecular symmetry-breaking and a transition from 
the normal mode to the local mode limit for the stretching vibrations due to rovibrational 
interaction. The application of the present theory with realistic molecular potentials to the 
H2Te, H2Se and H2S hydrides results in predictions of the bifurcation points very close to 
those calculated previously. However for the lighter H2O molecule we find that the bifurcation 
occurs at higher values of the total angular momentum than obtained in previous estimations. 
The present work shows it to be very unlikely that the bifurcation in H2O will lead to 
clustering of energy levels. This result is in agreement with recent variational calculations. 



1 INTRODUCTION 



Energy level clustering in molecules is very interesting phenomenon, which has been studied for 
about twenty years [|, |, | §. However the main interest was concentrated on highly symmet- 
rical molecules like tetrahydrides. The possibility of fourfold cluster formation in the upper part 
of rotational multiplets of symmetrical non-linear triatomic molecules has been predicted in the 
theoretical work by one of the authors and B.I. Zhilinskii |J. In that work the classical rovi- 
brational dynamics of an asymmetric top molecule was investigated in the rigid-bender model 
@. It has been shown that with the increasing of the angular momentum quantum number J 
the local precession around the axis with the smallest moment of inertia changes to delocalized 
precession around two equivalent axes at certain critical value J c . This phenomenon can be 
refered to as a bifurcation. It manifests itself in the softening of the precessional mode, i.e., 
in the decreasing of the spacing between the upper levels of the rotational multiplets when J 
passes through the critical value J c . The noticeable softening of the precessional modes in the 
ground and v<i vibrational states of H2O and H2S was interpreted as a possible bifurcation in 
this molecules || 0]. More precise numerical calculations based on the quantum Hamiltonian 
of the rigid-bender model were carried out soon ||, g]. They demonstrated the formation of 
fourfold energy clusters for J greater that J c . The microscopic theory of bifurcations developed 
in Ref. || allowed the estimation of the critical angular momentum J c which determines the 
region of cluster formation. It was found in Refs. [0, [l(]] that J c decreases with increasing mass 
of the central atom X in the H2X molecule. This fact has been used in the experimental study of 
"level clustering" in the rotational spectrum of H^Se [|ll|, 12 1. It was established that the energy 
levels in the vibrational ground state of H^Se form groups of four quasi-degenerate levels at the 
top of the J-multiplets when the value of the total angular momentum J exceeds the critical 
value J c . This phenomenon is very unusual for the rotational spectrum of a rigid asymmetric 
top. 

Soon afterwards the MORBID program has been successfully applied to the cluster analyses 
of the spectra of H^Se [13, 14], H2S 15] and H^Te |l6|]. The MORBID program is a variational 
approach allowing the calculation of energy levels of the total rovibrational Hamiltonian of a 
molecule in an isolated electronic state [17|. MORBID calculations based on fitted or ab initio 
molecular potentials showed good agreement with known experimental data and predicted the 
formation of fourfold clusters in highly excited rotational states of H2Se, H2S and H2TC The 
MORBID calculations were also in good agreement with pure classical, semiclassical and model 
quantum mechanical estimations of the critical angular momentum J c for these three molecules. 
However preliminary calculations [|l^] of the rotational levels up to J = 42 in the ground vi- 
brational state of water molecule using the MORBID program with optimized potentials [19 
did not show evident fourfold clusters. Recently Polyansky et al. have reported a new opti- 
mized water potential [2C[. Their calculations up to J = 35 for the ground vibrational state 
reproduce very well the available experimental energy levels but the authors concluded that the 
"rotational energy level structure in water is at least of a very different nature than the fourfold 
cluster structures observed for H2Se and calculated for H2S, H2Se and H2Te" [pQ|| - These two 
facts are in disaccord with all estimations of J c for the water molecule: J c = 27—28 in Ref. Q, 
J c = 26 in Ref. || and with the calculations in Ref. |2l]] , which predict quasi-degeneracy of the 
levels Jj j0 , Jj,i, Jj-x,i, and Jj-1,2 for J > 32. 

This extensive calculations show that the role of stretching vibrations in the bifurcation 
phenomenon is not clear enough. The rigid-bender model used in the early microscopic theory 
|| assumes the bond lengths to be fixed at their equilibrium values. This theory cannot also 
answer the question if a bifurcation exists in the rotational bands of excited vibrational states. 
This deficiency is remedied in the present paper, which uses the fully-dimensional description 
of rovibrational motion in a triatomic molecule, the stretching vibrations being included. The 
analysis of the classical equations of motion is based on the study of small harmonic oscillations 
near stationary states of the system. Four-dimensional oscillations can be separated into slow 
precessional and fast vibrational motion. The approach of the present work differs from the stan- 



dard adiabatic approximation by Born and Oppenheimer in that it considers the rovibrational 
motion to take place near the stationary states of a fast uniformly rotating molecule instead of 
having them take place near the equilibrium state of a non-rotating molecule. The difference 
has very important consequences. First, the precessional motion of an AB2 molecule around the 
axis with minimal moment of inertia becomes unstable at the bifurcation point J c . For higher 
angular momenta, the fast uniform rotation of a molecule takes place around two axes situated 
in the molecular plane between the axes of maximal and intermediate moments of inertia. Ow- 
ing to molecular symmetry, these two axes are equivalent. The delocalized precession around 
them may result in fourfold level clustering. This effect is the spectroscopic manifestation of 
the bifurcation. Another consequence is the asymmetrical deformation of the AB2 molecule by 
centrifugal forces, whose anisotropic action on the B-atoms causes one A — B bond length to be 
longer than the other. This symmetry breaking of the molecular configuration in its turn changes 
the vibrational dynamics in such a way that the normal stretching modes v\ and vj, transform to 
the local vibrations either of two A — B bonds. It is important to note that the transition from 
normal mode to local mode vibrations is caused in our case by rovibrational interaction and not 
by the anharmonicity of the vibrations as in customary local mode theory (see, for example, the 
review article (2^] ) . The developed theory allows to understand the relationship of local mode 
vibrations to a cluster formation in rotational spectra obtained by a model consideration in Ref. 




The classical method developed here has been applied to the hydrides H2O, H2S, EbSe, and 
H^Te to demonstrate that it produces consistency with purely quantum mechanical calculations 
for quantities like the critical angular momentum, the bond length distortions and the bending 
angle of the rotating molecule. The results obtained for the EbTe, H2Se, H2S molecules are very 
close to those obtained previously, but for the lighter H2O molecule the bifurcation is found to 
occur at J values significantly higher than those obtained in previous estimations. The most 
probable explanation of this is the neglect of the stretching vibrations in the earlier models. 
Furthermore, our analysis shows also that the stretching vibrations result in drastic reduction of 
the potential barrier, which separates the two regions of a delocalized precession. Thus it is very 
unlikely that the bifurcation in the vibrational ground state of H2O can result in an observable 
fourfold cluster structure. The clustering of levels cannot therefore be the ultimate criteria of 
the bifurcation. 



2 CLASSICAL EQUATIONS OF ROVIBRATIONAL MOTION 
FOR AN AB 2 MOLECULE 

We begin with the classical rovibrational Hamiltonian which can be derived following the method 
of Wilson, Decius, and Cross [g3]. Our treatment differs from theirs in that we do not use normal 
coordinates to have the possibility of considering the large displacements from the equilibrium 
configuration. With the standard choice of Euler angles |2^| as rotational coordinates, the 
classical kinetic rovibrational energy in a rotating coordinate system with the origin in the 
center of mass of molecule has the form 

T =\j2 kjUiWj + Ui 9i + 2 H m <xvli , ( 1 ) 

ij i ai 

where the indices i and j run through the axes x, y and z of this coordinate system, r a (x a ,ya, z a ) 
is the position vector of nucleus a with mass m a , v Q is its velocity, and g = Yl a m a( r a x v o)- 
Iij is a matrix element of the inertia tensor, and uJi are the projections of an angular velocity 
on the axes x, y and z: 

oj x = sinx — if sin cos X, 
u)y = 9 cosx + <£> sin sin x, 

uj z = tpcosO + x- (2) 



We can further introduce the total angular momentum vector J with projections 



Jx = Pe smx — Pip csc6>cosx + Px cot # cos 
Jy = Pe cosx + JV esc sin x — p x cot 9 sin x, 

Jz = P x , (3) 
which is related to the angular velocity projections as follows [24] 

J i - 9i = ^kjUj- (4) 



The inversion of Eq. (||) leads to 



u 3 =^2^ij( J i- 9i), (5) 



where {nij} is the matrix inverse to {Iij}- As a result the kinetic energy takes the form 

I I 

T = g J2 ( J i J i ~ 9i9j + 2 m * v li ■ ( 6 ) 

ij ai 

We consider a triatomic molecule B1A2B3 with nuclear masses m\ = 7773 = m and mi = M. 
To define the rotating (molecule fixed) axis system we can use one of the methods suggested by 
Sutcliffe and Tennyson ]25fl . Let the molecule be placed in the (xz) plane, i.e., yi = 0. Then 
we define the x-axis to be parallel to the bisector of the bond angle a = B1A2B3. The axis 
directions are chosen so that z-axis points from B3 to Bi, x-axis points from the centre of the 
mass to A2 and the coordinate system is right-handed. It is natural to introduce three internal 
coordinates: the A2-B1 distance T\ = q%, the angle a = qi and the A2-B3 distance 7-3 = q^. 
These internal coordinates are connected with the cartesian coordinates xyz as follows 

777x3 — {M + m)r\ a 

xi = — cos—, 

M + 2m 2 ' 

m(r± + rs) a 

X2 = — z~z COS—, 

M + 2m 2 ' 

mr\ — r^(M + m) a 

X3 = sm — , 

M + 2m 2 

(M + m)r\ + 7-3777, . a 
z\ = sm — , 

M + 2777 2 

m(r\ — 7-3) a 

22 = ~ M + 2m Sln 2' 

(M + 777)r3 + 7*1 777 a 

23 = wi^> sm o- 7 ) 

M + 2m 2 

With a molecule-fixed axis system defined in this way, the non-zero components of the inertia- 
tensor can be written as 

m(M + m) 2 2 . 

J-xz = — ttt r(?"i — 7-q sma, 

2(M + 2m) v 1 3J 

IXX = M + 2777 [M( ^ + ^ + m(ri + r3)2] Sln2 f ' 
IZZ = M + 2777 + r ^ + m(ri ~ r3)2] C ° S2 f ' 

m(M + 777) n 2 . 2m 2 

-f™ = — 7t (Ti + 7*0 ) — -— rir3Cosa, (8) 

yy M + 2m v 1 37 M + 2777 3 ' w 

and the elements of the inverse matrix are 



— n 



f^xx — 1 ' zz I (1 zz-^xx -^xz)i 
i^xz — ^xz/ij-xz Izz^xx)} 

[^zz — J-xx / {,1 zz^xx ^xz)' (*^) 

With the introduced internal coordinates the kinetic energy (^) can be transformed to the 
form 

T = JiJjVij + 2 ^Cvv'Qvqi/, (10) 

where i and j again assume the values x, y, z, and v and i/' run through 1, 2, 3. The c vv i 
coefficients have the form c vv > = a uu > — J2ij Gi v Gj v >^ij, but for a planar molecule they can be 
simplified to c vv < = a uv > — G yv G yu i \i yy . The matrix elements a uu > depend on internal coordinates 
as follows 

m(M + m) 

an = 033 - 



M + 2m 

m o , 2n 



022 = 4( M + 2m ) [( M + m )( r i +r 3 ) + 2mrir 3 cos a], 



2 

m 



° 12 = a21 = 2(M + 2m) r3Sma ' 



2 

m 



023 = a32= 2(M + 2m) riSma ' 



2 

m 



ais = aa^-^j^^ooaa, (11) 



and three G yv are written as 



2 

m 



^f 1 = ~o TT? r3sina ' 

2m + M 

m(m + M) 2 2 

Gy2 ~ 2(2m + M) (ri " r3) ' 

2 

G„3 = risina. (12) 

v 2m + M y J 

By introducing conjugate momenta p u = dT/dq u we can write the kinetic energy in the form 

T = JiJjfiij + -^2,b vv ,p v p v , , (13) 

ij vv' 

where {^w} = {c^^}" 1 assuming Det{c uu '} ^ 0. Finally the rovibrational Hamiltonian of the 
symmetrical AB2 molecule is expressed by 

H = -^Y^JiJjmj + ^Y^Kv'VvVv' +V(n,a,r 3 ), (14) 

ij vv' 

where V is the nuclear potential, which is symmetrical with respect to permutations of identical 
nuclei. 

The obtained Hamiltonian allows to write Hamilton's equations of motion 

Qv = ^K'vPv', (15) 
v' 

1 dfiij 1 db v > v " dV 



2 Y dq v 2^„ dq v dq v 

J% = ^^2 e ijkJjfJ-klJl- (17) 



The equation (|i~7| ) is derived by using the Poisson bracket {Ji,Jj} = eijkJk, where e is asym- 
metrical tensor and repeated indexes are summed. Equation ( |T7| ) is equivalent to the three 
equations 

Jx — [^xzJx x/J'Siz l^"yy)Jz\Jyi (1^) 
Jy — (l^xx H>zz)JxJz ^xz{.J z ^x)i (^) 
Jz — [ i^xzJ z \fJ"yy l^xx)Jx\Jy (20) 

They are not independent since J 2 is an integral of motion. 

3 PRECESSIONAL AND VIBRATIONAL MOTIONS OF AN 
ROTATING AB 2 MOLECULE 

Our analysis of the equations of motion is based on the study of rotation-vibrational motion 
near stationary states. In a stationary state, the time derivatives given in the left hand sides of 
Eqs. (lS-O) are equal to zero. Hence, in a stationary state Eq. ( |l5| ) takes the form 



J2 h v'vVv> = 0, (21) 

and it is easy to see that p u = 0, since the determinant of the {b uu /} matrix is non-zero as 
mentioned above. Now the right part of Eq. (^) can be written in the form 

= 0, (22) 

where the effective potential V e f[ is given by 

V e s = ~Y,Hi j JiJ j + V. (23) 

ij 

The three Eqs. (]22]) determine the configuration (i.e., the values of r s ±, a s , and r S 3) of a rotating 



molecule in a stationary state. The effective potential V e ff of Eq. (23) includes the centrifugal 
energy, which plays a central role in the dynamics of a rotating molecule. In particular, the 
configuration of a rotating molecule in a stationary state is different from that of non-rotating 
one. We will distinguish the corresponding values by indexes s (for the stationary state) and 
e (for the non-rotating molecule in its equilibrium configuration), respectively. Although the 
function [i xz is antisymmetrical under the permutation of the coordinates r\ and r% (see Eq. 
(|8|)), the effective potential (Eq. (|23"1)), the total Hamiltonian (Eq. (|l4])), and the corresponding 
equations of motion are invariant under the operations of the C2u(M) permutation-inversion 
group of the AB2 molecule (2(| (see also Ref. |p7j ). This group includes identity operator E, 
the P(13) permutation operator of identical nuclei, the operator of inversion E* and the product 
P(13)E* (see Table |. 

Equations (|l8|-pc|) together with the conservation law for J 2 define the axis of uniform 
rotation of the molecule, i.e., the total angular momentum vector J s in the molecular frame. 
Three projections of this vector are determined by four algebraic equations. Therefore, in a 
stationary state J s cannot have an arbitrary direction. There are two types of stationary states. 
In three axial stationary states Sf J s i = ± J, i = x,y, z, the molecule rotates uniformly around 
the z-axis. The molecular equilibrium configuration in the state Si is defined by the three 
equations 

2 dq v dq u 

with v = 1, 2, 3. The molecule has a symmetrical configuration with r s \ = r s ^ = r s in an axial 
stationary state Sj. The energy of this state is 



Ei = -it,;,;(r s , a.,, r.,) J 2 + V(r.„, a.,, r„). 



(25) 



There is only one plane stationary state S xz . In this state, the molecule rotates uniformly 
around one of the two equivalent axes situated in the (xz) plane symmetrically relative to the 
x-axis. They form angles (3 S and (tt — /3 S ) with the z-axis. The configuration of the molecule 
and the value of (3 S are defined by the equations 

sin p s + — sin 2p s + — — cos p s I J + — — = 0, f = 1,2,3 (26) 



2 V dg,, ' dq u dq u J dq u 

-xifJ-xx ~ Vzz) sin 2f3 s + fj, xz cos 2f3 s = 0. (27) 

It is interesting that in the state S xz , the molecule has an asymmetric configuration with r s \ ^ 
r s %. Of the two bond lengths, the one which initially forms the largest angle with the rotation axis 
will experience a larger elongation due to centrifugal forces, and for increasing J this bond length 
will tend to become perpendicular to the rotation axis. When we have obtained one solution of 
Eqs. (p6j-p7[), the remaining three solutions can be easily obtained by using symmetry operations 
of the C2v(M) group from Table |. The energy of the molecule in the S xz state is equal to 

E xz = ^J 2 (/i^sin 2 /3 s + /j, xz sin2(3 s + \i zz cos 2 (5 S ) + V(r sl ,a s ,r s3 ). (28) 

The analysis of stationary states shows that an axial state has higher symmetry in the C2«(M) 
group than the plane state. This means that transitions from axial to plane state or vice versa 
are accompanied by a C^-type bifurcation pp[ . Such a bifurcation has been considered for the 
first time in rotational spectra of symmetrical triatomic molecules in Ref. [0]. However, with the 
rigid-bender model with fixed bond lengths used in that work, it was not possible to show that 
the bifurcation results not only in the splitting of the stationary axis and consequently in level 
clustering but also in molecular symmetry-breaking. Due to this latter effect, the vibrational 
dynamics of the rotating molecule changes as we shall see below. Table || shows the centrifugal 
distortions Ar u = r su — r e of the bond lengths in the stationary state S xz . The values were 



obtained numerically for the molecules R^Se and H2S by using the potentials from Refs. [13] and 



[15], respectively. For comparison, the Ar u values obtained for the quantum cluster states in 
the upper part of the J-multiplets belonging to the ground vibrational state of these molecules 
are given in the same table. It is seen from Table [II] that classical and quantum estimates are 
in a good agreement. 

Since the derivative q^ z is positive in the stationary state S z , the bending angle a s decreases 
with increasing J according to Eq. (|24|). After a bifurcation point, the angle a s is stabilized 
at a constant value in the stationary state S xz . In the rigid-bender model this angle is inde- 
pendent of J for J > J c H (see Fig. |l|). For H2Se, the rigid-bender and stretching-bender 
models give comparable results, but for H2O the predictions of the stretching-bender model 
deviate significantly from those of the rigid-bender (Fig. |l|). In the stretching-bender model, 
the bending angle value for H2O continues to decrease significantly when J is above the critical 
value. Hence no stabilization at a constant value takes place. This effect can be explained as a 
competition between centrifugal force causing the hydrogen nuclei to bring together and their 
mutual repulsion. 

In order to follow the change of molecule rotational regimes as J increases, we must inves- 
tigate the stability of the stationary states. To do this we consider the linearized set of Eqs. 
(|l5T-ji"7|) for small displacements of the internal coordinates Q u = q v — q su and the angular mo- 
mentum projections J- = Ji — J s i from their stationary values. Let us begin with the axial 
stationary state S y . The linearized equations have the form 

W n dV cff \ 
\ dq u dq u i J 

Jx J {/J'xx H'yy){P'ZZ ^yy)'Jx — 0' (29) 

The coefficients of these equations are assumed to be taken in the stationary state, i.e., for 



state. Equations (29) describe two independent harmonic motions: vibrational and precessional. 
The latter motion is characterized by the precession frequency 



— J[{^xx Hyy ) (.Hzz Ms/j/)] ^ > 



(30) 



which has the same form as for a rigid top ||2 
differ from that of non-rotating molecule except in that the frequencies ui, ui2, and L03 of three 



The vibrational motion in the state S y does not 



normal modes vi, vi, and 1/3 [24] are determined by the effective potential and the equilibrium 
configuration of rotating molecule. As a result they are shifted relative to those of non-rotating 
molecule (see below). 

The linearized equations for the stationary state S z 



sr ( a d 2 V eS \ d[x xz 

y. c vu 'Q u i + - — - — Q v > + J— — J x = 0, 

<-f \ dq v dq v i / dq v 



Jx ~\~ J {p>zz Hyy) i {Hzz Hxx)Jx J ^ ] o Qv \ — ^' 

I v 0( lv J 



(31) 



describe coupled vibrational and precessional motions. It is not difficult to see that the coupling 
originates in the antisymmetrical variable Q a = (Q\ — Q^)/\/2. Thus Eqs. (31) split into two 
independent parts for the variables Q a , J x and Q s = (Q± + Q^)/V2, Q2- The equations for Q s 
and Q2 are independent of the precessional motion and describe the v\ and V2 vibrational modes 
modified by molecular rotation. We can separate the precessional motion from the antisymmetric 
vibration by using the inequality W3 3> Vt z (adiabatic approximation), which is correct up to the 
critical point J c (see below). In this approximation, the frequency of molecular precession around 
the z-axis is equal to 



j 



{fJ'zz Hyy) \ Hzz Hxx J 



( dn a 



V dr a 



< d 2 V eS 
dr 2 



1/2 



(32) 



where according to Eqs. (||) and (J9|) 



Hi 



Hxx — 



2(M 



m 



mMr 2 sin a s 



m 



M + m 



cos a. 



(33) 



For the hydrides considered, the equilibrium angle a e is larger than 90°. So the precession 
frequency of Eq. ( |32|) is real and the stationary state S z is stable for small values of J. As we 
saw above, the stationary angle a s decreases with increasing J due to a centrifugal distortion. 
So [i zz approaches [i xx and the precession frequency Q z vanishes when 



Hz 



1 j2 

Hxx 1 <J c 



QHxz 

dr a 



' d 2 V eS 
dr 2 



0. 



(34) 



Equations ( |34| ) and ( p4j ) (for i = z) define the critical angular momentum J c . In the rigid-bender 
model, this point is determined by the bending angle 



Or 



arccos 



m 



M + m 



(35) 



for which I xx = I zz . The last term in Eq. ( |34j ) takes into account the non-rigidity of the bond 
lengths r\ and and becomes important when \i zz — fi xx is close to zero. Since the correction 



term is positive, Eq. (34) gives higher values of J c than the rigid-bender model. Clearly, the 
more rigid the antisymmetrical stretching mode 1/3, the closer is J c to the value obtained in Ref. 
in the rigid-bender approximation. Our numerical calculations of critical angular momentum 
for the set of the H2X hydrides with the potentials from Refs. (T 
corrections to the J c values predicted previously. We obtained J, 



15, [l(| give mainly slight 
9.3 for H 2 Te, J c =12.5 for 



,A 7 _1 Q O f^v XT C 



on t~\ d~r 



,fS 1; 11/1 



It is seen that for lighter molecules, larger discrepancies are obtained. An exception is the water 
molecule for which the discrepancy of the J c value is very significant. Using the potential from 
Ref. J c =26.5 was obtained with the "frozen" bond lengths (which is close to all previous 
estimations) and J c =35.2 in the stretching-bender model. 

For J > J c , the stationary state S z corresponds to an unstable saddle-type point. It is the 
top of a separation barrier between two symmetrical states S xz which appear because of the 
bifurcation. The equations of a small amplitude motion near this state have the form 




n i d2y ^ n 

Oq v aq v i 
— Hyy + fJ-xz tan f3. 
where the "coupling constant" 

B v = 




0. 



J' x - J cos (3 S J2B V Q V ) = 0, 



1 % 
2 



t'z 



dq u 



■ sin 2fj s + — — cos 2p s 



(36) 



(37) 



mixes all four motions: three vibrations and a precession. The last one can be decoupled from 
the vibrations in the adiabatic approximation for moderate values of J. In this approximation, 
the frequency of the precession around the S xz state is 



j 



(fJ-zz - Hyy + fJ-xz tan (5 S ) 



J 



cos 2j3 s 



+ x E(-ir Ku>b v b v 



1/2 



(38) 



where the determinant A = Dct 



\dq v dq v i J 



and A, 



are the minores of this determinant. 



All the values in Eq. (35) are assumed to be taken in the stationary point determined by Eqs. 
( |2"6| ) and (27). It is not difficult to show that the precession frequency Sl xz vanishes in the critical 
point J c . 

The small amplitude motion near the state S x is described by Eqs. ( |3l|) when index x is 
replaced by z. As a result we get in the adiabatic approximation the following expression for 
the precession frequency around the x-axis 



a, 



j 



d^xz 
dr a 



' d 2 V cS 
dr 2 a 



1/2 



(39) 



where all the values are taken in the stationary state S x . The last term in the braces is a small 
correction. Consequently, the stability of the S x state is mainly defined by the value given 
in Eq. ( |33"[ ) which in turn depends on the angle a s . Since this angle is larger than 90° for the 
molecules considered and increases with increasing J because the derivative -ffi 2 is negative, the 



precession frequency Q x is imaginary for reasonable values of J. Thus the state S x corresponds 



to an unstable saddle-type point. The J-dependence of its energy is defined by Eq. (25). 

Equations ( pq ) and ( p8| ) (when the transition S z — > S xz occurs) provide a good estimation 
of the maximal rotation energy in J-multiplet. The difference between our estimations and 
experimental data or variational calculations for H^Te, E^Se, H2S and H2O molecules does not 
exceed 10% . This is also true for our estimations of the minimal rotational energy in the 
J-multiplet using Eq. ( |25| ) for the state S y . The energy of the unstable stationary state S x 
corresponds to the top of the barrier responsible for the inversion doublets (-KT-doublets) in the 
upper and lower parts of the J-multiplets p. Another barrier appears because of the bifurcation 
for angular momentum at J > J c . It separates two degenerate maxima of the stationary state 
S xz . The barrier is responsible for the twofold clustering of the inversion doublets. Its height 
is equal to the difference between the energies of the states S z and S xz . The higher the barrier 
is, the more apparent is the clustering. In the EbTe, H^Se and H2S molecules, the height of the 
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with potentials from Ref. show that the height of the energy barrier increases very slowly 
with increasing J and its value is only 15 cm -1 at J =40 (see Fig. ^). This value is less than the 
precession frequency in the S xz state so that such a barrier cannot result in fourfold clustering 
of levels. For higher J values the barrier increases further to its maximal value 150 cm' 1 (which 
is about the precession frequency at this point) and then begins to decrease. In contrary, the 
barrier rapidly increases with increasing J for the rigid-bender model as it is seen from Fig. |2[ 
where the dependence of the barrier height on angular momentum J is plotted. It is obvious 
that water represents a special case in the series of H2X hydrides. Despite the bifurcation in 
the ground vibrational state of this molecule, no fourfold level clustering can be observed in the 
upper part of J-multiplets since the barrier is too small. This fact gives us a physical explanation 
to the recent numerical calculations of water rotational spectrum with highly accurate optimized 
potentials [|l^, 2C], which failed to show any fourfold clusters. 

The C21, bifurcation has to exist also in rotational energy structure of excited vibrational 
states so far as the adiabatic approximation u 3> Vt is valid and allows to separate the precession 
from the vibrational motion. As to the energy clustering effect, we have seen above that this 
problem requires separate and more careful consideration. Preliminary analysis of experimental 
data showed the tendency of levels to cluster in V2 vibrational states of the H2O and H2S 
molecules |5|, 0. The cluster phenomenon has been found in variational quantum calculations 



for the excited vibrational states of H^Se [jllj, 14], H2S ||15[ and H^Te |T6| and very recently 



observed experimentally in the v\/v% state of E^Se [p9[|. 

Consider now in detail the molecular vibrational motion near the stationary states S z for 
J < J c and S xz for J > J c . Four normal modes of precession- vibration motion are described 
by four coupled Eqs. ([H]) and (36) correspondingly. All these equations are not correct near 
the critical point J c where the rovibrational motion is strongly non-linear. It was shown above 
that the equations for the axial state S z (and those for the state S x ) can be separated into two 
independent pairs when symmetrized stretching coordinates are introduced. When J= the 
precessional motion is absent and the vibrational motion corresponds to the well-known picture 
of the normal vibrational modes v\ , v 2 and 1^3 p4| . This limit allows us to refer to the standard 
notations for a rotating molecule as well, i.e., the frequency u)\ corresponds to the symmetrical 
stretching mode u±, the frequency L02 corresponds to the bending mode f2, and the frequency 
U3 corresponds to the antisymmetrical stretching mode V3. The difference from a non-rotating 
molecule consists in the effective potential for vibrations and in the equilibrium configuration, 
which is symmetrical but differs from that of non-rotating molecule. Besides, the precessional 
mode with frequency £l z appears and the antisymmetrical stretching vibration mixes with the 
precessional motion because of a rovibrational interaction. However this mixing is very weak. 
This can be seen from Table III where the normalized amplitudes of four coordinates Q\, Q2, 
Q3, and Q4 = J x /J involved in a normal mode oscillations are given for the H^Se molecule at 
J= 10. It is seen from the table that the mixing of different types of motions is small and the 
character of the vibrations is close to the standard normal mode picture p4]j . 

The character of the oscillations changes drastically after the critical point J c . It is impos- 
sible to split the four equations ( j36[) describing the small amplitude motion near the stationary 
state S xz into two independent pairs because the molecular configuration is asymmetric and the 
uniform rotation of molecule around the axis in the (xz) plane mixes the symmetrical and anti- 
symmetrical stretching vibrations. As a result, the standard picture of normal mode vibrations, 
which is correct for the S z state, breaks down. This effect is illustrated in Table IV which shows 
normalized amplitudes of oscillations around the asymmetrical configuration of the state S xz 
at J=16 {(5 S = 26.08°). It shows also that the precessional motion is mixed mostly with the 
bending vibration. It is worth noting that the vibrations in the mode with the frequency u\ are 
localized predominantly on the bond r% and vibrations in the mode with the frequency L03 are 
localized on the bond r%. This feature is inherent in local mode vibrations p2| . However usually 
the local mode motion appears due to the strong anharmonicity as the stretching vibrations are 
increasingly excited. In our case the normal to local mode transition is a result of the symmetry 
breaking, which leads to the different force constants for two bonds and favors the local mode 



picture in the competition between kinetic and potential coupling between the local bond oscil- 
lators. Accordingly the local character of the stretching motion strengthens with increasing J. 
This follows from Table [V| which shows normalized amplitudes of oscillations around the state 
S xz at J=40 (p s = 42.26°). 

The classical calculations of the present work have shown that in non-linear AB2 molecules, 
a normal to local mode transition can be brought about by rotational excitation. Similar results 
have been obtained in quantum numerical calculations with the MORBID program [14]. The 
classical analysis of precession- vibrational motion of the present work gives a physical picture 
of the phenomenon. That is, the bifurcation in a ground vibrational state breaks the molecular 
symmetry and changes both rotational and vibrational motions of the non-linear AB2 molecule 
resulting in the fourfold clusters and in the normal to local mode transition. The physical 
reasons for these phenomena lie in the molecular symmetry and in the anisotropic centrifugal 
force. Therefore we should expect similar or even more interesting phenomena in the fast rotating 
molecules with higher symmetry. 
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FIGURE CAPTIONS 



Figure 1: J-dependence of the bending angle in the stationary states with maximal rotational 
energy (S z if J < J c and S xz if J > J c ) for the H^Se and H2O molecules. The angle is displayed 
relative to the equilibrium value a e for non-rotating molecule. The solid lines correspond to the 
stretching-bender model with the critical angular momentum J c =12.5 (H^Se) and 35.2 (H2O); 
the dashed lines represent the rigid-bender model (J c =11.4 for the H^Se and 26.5 for the H2O 
molecules) . 

Figure 2: Value of the barrier which separates two equivalent maxima of stationary state S xz 
for water molecule calculated in the stretching-bender model (solid line) and the rigid-bender 
model (dashed line). 



TABLE CAPTIONS 



Table I: Action of C2«(M) group operators on the internal coordinates and the projections of 
the total angular momentum J. 

Table II: Comparison of the molecular centrifugal distortions Ar u = r su — r ev for the stationary 
state S xz at J =20 calculated in quantum and classical methods. The classical calculations cor- 
respond to molecular rotation around the stationary axis which is approximately perpendicular 
to the bond r±. Indexes s and e stand for the molecular equilibrium configurations of rotating 
and non-rotating molecule correspondingly 

Table III: Normalized amplitudes of coordinate oscillations in the normal modes of precession- 
vibrational motion (Q4 = J x /J) around stationary state S z for H^Se molecule at J=10. 

Table IV: Normalized amplitudes of coordinate oscillations in the normal modes of precession- 
vibrational motion (Q4 = J' x / J sz ) near the equilibrium configuration of the H^Se molecule 
uniformly rotating in stationary state S xz at J = 16 around the axis which is approximately 
perpendicular to the r\ bond (f3 s = 26.08°). 

Table V: Normalized amplitudes of coordinate oscillations in the normal modes of precession- 
vibrational motion (Q4 = J' x /J sz ) near the equilibrium configuration of the E^Se molecule 
uniformly rotating in stationary state S xz at J=40 around the axis which is approximately 
perpendicular to the r\ bond (/3 S = 42.26°). 
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Table 0. 



molecule 


quantum estimation 
Ari/ A Ar 3 / A 


classical estimation 
Ari/ A Ar 3 / A 


H 2 Se 
H 2 S 


0.023 a 0.004 a 
0.024 6 0.009 6 


0.025 0.002 
0.022 0.010 



a — taken from Ref. Q . 
b — taken from Ref. gj|. 
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III. 



Normal 
frequency /cm -1 


Normalized amplitudes 


(n - r s x)/k 


(a — o s )/rad. 


(T3 - r s3 )/k 


J x /J 


u t = 2426.4 


l.t 


-0.032 


1. 


0. 


u 2 = 1082.3 


-0.014 


l.t 


-0.014 


0. 


uj 3 = 2437.4 


l.t 


0. 


-1. 


-0.004 


tt z = 15.7 


0.004 


0. 


-0.004 


l.t 



Normalization condition Q n =l. 



Table IV. 



Normal 
frequency /cm~ 1 


Normalized amplitudes 


(n - r s i)/k 


(a — a s )/rad. 


(T3 - r s3 )/k 


yj x J sx) / J sz 


ui= 2390.9 


l.t 


0.018 


0.095 


0.003 


L0 2 = 1110.1 


-0.015 


l.t 


-0.013 


0.020 


lo 3 = 2436.0 


-0.095 


0.012 


l.t 


-0.003 


Q xz = 39.5 


0.006 


0.081 


-0.008 


l.t 



Normalization condition Q n =l. 



Table 0. 



Normal 
frequency /cm~ 1 


Normalized amplitudes 


(n - r s {)/k 


(a — a s )/rad. 


(T3 - r S 3)/A 


\J X J sx) / J sz 


u x = 2095.8 


l.t 


0.030 


0.015 


0.004 


lo 2 = 1271.9 


-0.022 


l.t 


-0.009 


0.093 


lo 3 = 2438.0 


-0.015 


0.003 


l.t 


-0.002 


tt xz = 211.7 


0.000 


0.429 


-0.009 


l.t 



Normalization condition Q n =l. 



